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Abstract. We propose to contruct a temporary wave on the surface of the ocean, as a particular 
sohition of the Saint- Venant equations with a source term involving the friction, whose shape is 
expected to mimic a rogue wave 

The phenomenon of Rogue Waves, or Freak Waves, is a transitory phenomenon which appears 
in the open ocean, under the shape of a gigantic and devastator wave, because it accumulates 
a significant quantity of energy. Many articles are devoted to them in [6]. We propose here an 
hydraulic model, making use of the Saint- Venant equations. The model is limited to only one 
dimension of space, in the direction of the propagation, since the side movements near the wave are 
uniform. The validity of this model requires a large wavelength, particularly in deep ocean. The 
wave itself has a rather short wavelength, but the whole phenomenon involves two waves of large 
wavelength, and relatively low amplitude. 

We denote by H the mean depth of the ocean, and by the sonic velocity of waves inside 
the water (we have Cg = 1647 ms~^). The wavelength A of the Saint- Venant waves must satisfy a 
condition of the form 



0.1 A > 2 — , 

Cs 

where g is the gravity constant and the number of sonic interactions (back-and-forth) between 
the bottom and the surface of the ocean. This condition means that along a horizontal distance of 
a wavelength A, there are at least such sonic interactions . The use of the Saint- Venant model 
is as more appropriate as is great. We usually require A^ > 25, which implies for example a 
wavelength greater than 21400 m for an ocean depth of 3700m. We suppose here, for simplicity, 
that the propagation goes from West {x < 0) to East {x > 0). 

1. The INITIAL CONFIGURATION OF THE WAVE 

We denote by g = q{x, t) the ocean depth at a point x at a time t, and m = Tn(x, t) the corresponding 
flux. We have m = q u, where u is the water velocity. The relative velocity of the waves of the 
ocean surface is given by 

c = y/gq (= c{q)) . 

We denote by > the friction coefficient of Strickler and we consider as a simplification purpose 
that the bottom is fiat. In this configuration, the Saint- Venant model reads 

(1.1) + m, = , 

(1.2) rrit + 2u rrix + (c^ — u^) q^,. + k \u\u = . 
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We shall use the representation of waves proposed in [3] when a source term (here k\u\u) 
is present. The different states of the same wave are described by a segment of a straight line 

(1.3) m = Aq - B , 

in the phase plane, that is the plane (g, m) here. The parameter A is a constant which corresponds 
to the wave velocity, of the form A = Uref — c{qref) or ^4 = Uref + c{qref) depending if the wave 
is travelling eastwards (sign +) or westwards (sign — ) and for a given reference state Mref = 
{qref, TTTref), with TTiref = qrefUref , of course. The parameter B is also a constant and is determined 
by the reference state since B = rriref + A qref - We denote by Mq the state of the ocean far on 
the west side, and by the state of the ocean far on the east side. In both cases, the distances 
are supposed to be larger than the reference wavelength A proposed in (10. II) . which allows to use 
the Saint- Venant model. We suppose that these states also corresponds to reference velocity equal 
to zero; this assumption is linked to the hypothesis of a flat bottom of the ocean. That way Mq 
represents a state Mq = (go, 0) and represents a state = (g^,, 0). We suppose 

go > • 

A difference go — g* of a few decimeters is enough even for a wide depth of the ocean. The profile of 
the initial state is made of two branches which meet in a state P = {qp, nip) for example at x = 0, 
which corresponds to the choice of the origin. From the state P to the state M^, that is on the 
East side, the profile is decreasing and referred to the state M^. The corresponding states are so 
situated on the straight line 

m = {q- gj , 

with = c{q^) = ^/gq^ . As a matter of fact, we have A = and B = c^q^. for this part of the wave 
profile. The explicit formulation of the profile is obtained, following [Ij, by inverting the profile 
relation 

(1.4) i)E{q) = X , 

where the index E stands for East. The inverse profile iPe{q) is determined by integrating 

^^(^^ = k\Aq-B\{Aq-B) ' ^^^^^^ = ° " 

We need to have g* < go < gp , in order to ensure an increasing profile on the West side, then a 
decreasing one on the East side. We get that way a positive flux: Aq — B = c^{q — q^) > 0, and 
by writing ^ = — we get 



thus 



i-e i+i+e 1 



k{l-if k{l-i) k \ 



M<1) = In (^^) +h<lP-<l) + 7^ ill - ■ 

k \qp-q*J k 2kq^ 

By inverting this function ipEiq) and using (ll.4p we get g as a decreasing function of x which is 
equal to qp when x = . For the West side of the profile, with the index W, the reference state 
Mref = {Qref, "rriref) uiust correspoud to a depth satisfying 

Qref > qp (> go) , 
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in order to ensure an increasing profile. The reference velocity associated witli tlie state M^ef tliat 
is 

Aef = + Cref , Wtth Cref = C[qref) = y/gQref , 

Qref 

is also the velocity of the West side profile of the wave and is positive. The straight line representing 
this West profile in the phase plane has the form 

m = rriref + Aref (q-qref) , 
and passes by the state Mq = {qo,0). Hence 

'rriref = Aref {Qref — Qo) ; Bref = Cref Qref ■ 

The invert profile is described by a function i/jw satisfying 

/' / \ ref WefQref i f \ n 

k [Arefq - Bref) 

with qo < qp < Qref ■ The initial West side profile of the wave is then obtained by inverting, for 

any x < 0, 

(1.5) i^wiq) = X . 

We set ^ 

^ = , ^0 = — — 5 Fref = — — , ( Froude number ) , 

Qref qref Cref Qref 

to obtain 



Mq) = f-C-2eo + ^#^ + ^^) , 

k{Fref + i)' V e-eo (e-eo)v 

Iref 

that is 

Qp-Q^ ,c.Q^ I \ , 9o / . /?-?o\ , {Qlef - Qo) i.Q - Qp) 



^3 



iI)w{q) = K 
with 



+ 2— {qp - q) + — (go - ^ref) In ( — — — ) + 

Qref Qref \QP ~ QO / 



'2 Qref Qref Qref \QP ~ QO J Qref (Q - Qo) {QP - Qo) 

1 



K 



k {Fref + 1)' 

The whole initial profile is now given by (ll.Sp for x < and by (11. 4p for x > 0. It corresponds to an 
increasing function q{x, 0) for x < and a decreasing one for x > 0, which is continuous at x = 
where its value is g = gp. 

2. The propagation of the wave 



The wave profile is expected to propagate eastwards, with the respective velocities Aref and 
which are different for each part West or East. Since the west profile moves lightly faster, the crest 
of the profile will move up, at the junction of the two parts. The left side of this crest corresponds 
to the continuity of the West profile, extropolated for depth values q going from qp to the maximal 
value qref- 

The right part of the crest corresponds to a discontinuity, that is a shock wave, whose location 
is imposed by the mass conservation. At any time t the water contained in the bump under the crest 
comes from the column of water of length {Aref — c*) t shifted since the initial time. We denote 
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Figure 1 . Propagation of a Rogue Wave 

by qi{t) and qr{t) the water depths on the left side (index I) and on the right side (index r) of the 
discontinuity, and by Xo{t) its position. We always have 

q* < qr(t) < qp < qi{t) < qref ■ 

As the West profile is moving with the constant velocity Aref, we get it by simply inverting for any 
time t, the relation 

for X < xo{t), with the function ifjw defined above. By the same way, the East profile is obtained 
by inverting the relation 

i^E^q) = X - cj 

for X > Xo{t). with the function ■i/'s defined above. For a given time t, the depths qit) and qr{t), and 
the shock position xo{t) are linked by three conditions, entailing three equations. The first equation 
says that qi{t) is the value of the West profile at x = xo{t), that is 

i^wiqiii)) = Xo{t) - Aref t ■ 

The second equation says that qr(t) is the value of the East profile at a; = Xo(t), that is 

^E{qr{t)) = Xo{t) - C^t . 

The third equation is given by the compatibility relation of Rankine-Hugoniot 



/oiN / u\ u\\ I qr{t) + qi{t) , ^re/ ^0 c* 

(2.1) (?,(*) - + - ^ = \ef - . 

which ensures the mass conservation. A dichotomy method running on the parameter Xo{t) allows 
the simultaneous determination of these three parameters. For practical purposes it is however 
easier to determine Xo{t) by checking directly the mass conservation. Let us consider two points Xi 
and X2 such that 

Xi + Aref t < Xo{t) < X2 + t . 

If it is not the case in a first choice, one can increase X2 or decrease Xi sufficiently. We denote by 
Mo the mass of water laying between the two points xi and X2 at the initial time: 
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pX2 rO pX2 

= / 9(2^,0) dx = Qwix.O) dx + qE{x,0) dx , 

J xi J x\ Jo 

where qw and qE are the respectives depths of the West and East profiles. The same mass Mq has 
to be found at any time t between the two trajectories of equations 

(2.2) x[it) = A^^f Iref^ref ^ ^ ^^^^^ ^ ^ 

qw (xi{t) - A^^ftj 

and 

(2.3) x'^it) = c, (1 — ^ -) , X2(0) = X2 , 

which means that the relation M{t) = Mq always occurs, with 

M{t) = / q{x,t) dx . 

Jxi{t) 

We remark that the mass M{t) may be explicitely calculated from the relation 

ritit) rq-iit) 
M{t) = q i/j'wiq) dq + q i/j'^iq) dq , 

where qi{t) = qw{xi{t)) , q2{t) = g£;(x2(t)) . This calculation involves only primitives of rationnal 
fractions. More simply, using the computing files giving qw on the West side of a point xq and qE 
on the East side of this point xq, one can compute 

pxo rx2{t) 
F{xq) = qw dx + qE dx , 

J Xi(t) J Xq 



and next determine xq = XQ{t) such that 

F{xo) = Mo 
by noticing that F is an increasing function. 



This last process has been used to compute the results shown on the figure above. For this 
example, we have taken q^ = 3700 m and qo = 3700.2 m, then q^ef at its maximal value, that is 
qref = 3731.6737 m here. This choice leads to the value of qp, which is qp = 3715.8087 m here. The 
friction coefficient was taken equal to = 0.45 . We have to notice that the value of the friction 
coefficient is strongly linked to the wavelength of the profiles since the smallest frictions give the 
larger wavelengths. For terrestrial hydraulic flows (rivers or estuaries for example) the values of 
the friction coefficient is far smaller, of the order of 10~'^ but corresponds to far less deep flows. It 
seems natural to consider that for larger depths this parameter has to be upgraded. At the time 
t = 1000, the shock amplitude qi — qr reaches the value 11 m, between the two cells adjascent to 
the computed position of xo{t) . The relative error on the mass is of order of 10~^ (mainly due to 
the meshsize), and the trajectories (12. 2p and (12. 3p were approached by a simple trapezoid formula, 
since the depth has a small variation along the trajectories Xi{t) and 2:2 (t) when they are chosen 
relatively far from the shock. 
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Figure 2. Propagation of a bigger wave 



By upgrading lightly the value of qo, one get more important shock amplitudes. For the next 
example, with qo = 3700.8 m, we obtain the values qref = 3763.8773, and qp = 3731.8248. The 
computed shock amplitude qi — qr is equal to 33.787 m. The wave crest culminates at more than 
50 meters above the sea level. 

3. Graphic interpretation 
We propose a graphic interpretation in the phase plane (g, m) using the picture below. 




Figure 3. The representation in the phase plane 



The states Af^, and Mq are represented on the g— axis. The line passing through corresponds 
to the East wave, of equation m = c*(q' — g*). The line passing through Mq corresponds to the 
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West wave, of equation m = A^efiq — Qo)- These two lines meet at the point P. The shock wave is 
represented by the segment MiMr, with Mi shding along the West wave line, from P to M^ef and 
sliding along the East wave line, from P to M^,. This segment MiMr is not exactly a straight one, 
and its equation is obtained from the Rankine Hugoniot condition (12.11) . The maximal amplitude 
is reached when the first event between Mi reaching Mref or Mr reaching M* occurs. After that 
time the wave is expected to collapse. The scaling of the picture has been strongly modified for 
readability, since in the real example all those lines are so close to one another that it is impossible 
to perceive any difference. 



4. Conclusion 

This study shows that the outbreak of these waves is due to a differential in the pressure field, 
resulting from two waves with different profiles, with close but different velocities and sufficiently 
large wavelengths. The water is pushed up from this pressure effect together with a friction effect. 
The later behaviour is not studied here: one expects that, when qi{t) reaches the la value Qref (which 
occurs in a finite time), a backwards wave will appear, with a negative velocity, and will provoke 
strong perturbations on the West profile near the wave crest, which will collapse soon. However 
this backwards wave will have a shorter wavelength, probably too small to fit up with the use of 
the Saint- Venant model, as described by the condition (10.11) . 

Another remaining work is to fit out the different parameters in order to be in accordance with 
real world observations. We also outline the strong sensitivity of the difference q^ — q^ on the hight of 
the wave crest, and then on the shock amplitude. A little more important difference between q^ and 
causes noticeably more important shock amplitudes. We have only proposed empirical choices 
of the parameters, in order to get realistic results showing that this study may be a suitable way to 
understand the shaping of Rogue Waves. We emphasize the friction plays here a fundamental role 
since it allows a linear behaviuor of the West and the East waves. 

Another idea to retain is a new example of the application of the notion of Source Waves after 
the Roll Waves in channels and rivers, the Tidal Bore waves in the estuaries, the surf waves on the 
shore near the beach, the hurricanes and tsunamis (see [5], [1]). 

The authors thank Michel Olagnon from Ifremer-Brest for some useful answers by e-mail. Read- 
ing for example p] or [3] is also very instructive for the description of the phenomenon of Rogue 
Waves and starting some bibliography research. It appears that a discussion on either the linear 
character or the non linear character of the waves is spreading. In this study we put together 
both characters since linearity is introduced through the source term, that is the friction term, 
with different parameters on the two sides of the wave, and the nonlinearity effect is present in the 
Shockwave, linking the two sides of the wave. 
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